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Abstract. Continuum driving is an effective method to drive a strong stellar wind. It is governed by two limits: the Eddington 
limit and the photon-tiring limit. A star must exceed the effective Eddington limit for continuum driving to overcome the 
stellar gravity. The photon-tiring limit places an upper limit on the mass loss rate that can be driven to infinity, given the 
energy available in the radiation field of the star 

Since continuum driving does not require the presence of metals in the stellar atmosphere it is particularly suited to 
removing mass from low- and zero-metallicity stars and can play a crucial part in their evolution. 

Using a porosity length formalism we compute numerical simulations of super-Eddington, continuum driven winds to 
explore their behaviour for stars both below and above the photon-tiring limit. We find that below the photon tiring limit, 
continuum driving can produce a large, steady mass loss rate at velocities on the order of the escape velocity. If the star 
exceeds the photon-tiring limit, a steady solution is no longer possible. While the effective mass loss rate is still very large, 
the wind velocity is much smaller 
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INTRODUCTION 

Continuum driven stellar winds are comparatively rare. 
The star has to exceed the effective Eddington limit for 
continuum driving to be effective at all, with the effective 
Eddington limit defined as the point where gravity is 
balanced by the continuum force including both electron 
scattering and other continuum processes (i.e. bound-free 
and free-free). Since most stars stay below this limit their 
winds are driven by line absorption of radiation rather 
than continuum scattering. 

However, continuum driving has two characteristics 
that may yet make it of great importance for the evolution 
of zero- and low metallicity stars: 1) Unlike line driving, 
continuum driving does not depend on the presence of 
metals. Rather it depends on the presence of free elec- 
trons as well as bound-free interactions, both of which 
can take place in any gas, irrespective of its composition. 
2) Continuum driving can propel a potentially far greater 
mass loss rate than line driving. Whereas line driving is 
limited to about lO^'* — 10^^ M© yr^^ because of line 
saturation, continuum driving can easily drive a hundred 
times more mass. In fact, the only limit to the mass loss 
rate of a continuum driven wind is the total amount of en- 
ergy that is available in the radiation field of the star. For 
line driven winds, this upper limit is irrelevant, since the 
mechanical luminosity of such winds is always much less 
than the radiative luminosity of the star. The high mass 
loss rate makes continuum driving ideal for explaining 



massive outbursts like the one rj Carinae experienced 
between 1840-1860 (Davidson & Humphreys 1 1]; Smith 
|7]). 

The main problem in producing a steady, continuum 
driven wind lies in the fact that the ratio of radiative 
flux and gravitational acceleration tends to be constant 
throughout the star. (Both scale with r^^.) Therefore, if 
the star is super-Eddington at the surface, it will be super- 
Eddington throughout its envelope, which leads to large 
scale instabilities as the envelope becomes gravitation- 
ally unbound. 

A potential solution to this problem is provided by 
the "porosity" of the material. This decouples radiation 
from matter as photons escape through the low density 
medium, reducing the interaction with the main body of 
the gas, which is tied up in optically thick clumps (Shaviv 

Owocki et al. f^, henceforth OGS, worked out analyt- 
ically what a porosity modified, continuum driven wind 
would look like as long as the star does not exceed the 
photon-tiring limit. Here we present numerical models of 
similar winds, also exploring what happens when a star 
does in fact exceed the photon tiring limit, which can not 
be described analytically. 



ANALYTICAL AND NUMERICAL 
BACKGROUND 

For our simulations we assume that the stellar atmo- 
sphere consists of a porous material, where optically 
thick clumps are surrounded by a low density, opti- 
cally thin gas. This will cause the photons to escape 
through the optically thin medium, moderating the in- 
teraction with the gas in the clumps. As a result, the 
coupling between radiation and matter is reduced. This 
effect decreases as the clumps themselves become opti- 
cally thin, making the interaction between radiation and 
matter density-dependent. As a result, the full force of 
the radiation field only comes to bear on the gas in the 
low density outer layers of the stellar atmosphere. 



Analytical solution 

The acceleration produced by radiation force on mat- 
ter through continuum interaction is given by: 
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with K the continuum opacity, F the radiative flux and 
c the lightspeed. If the medium is porous, the standard 
opacity has to be replaced by a new, effective opac- 
ity, which scales with the optical depth of the individ- 
ual clumps. This in turn depends on the density of the 
medium. Rather than giving all clumps the same optical 
depth, we assume a power-law distribution, which leads 
to: 
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with Up the power-law index and Tq = p /po, with po = 
1 / (/ifc) the density at which the largest clump has optical 
depth unity. This depends on the opacity and the porosity 
length: h (OGS), which is related to the size of the 
clumps relative to the distance between clumps. 

Since the total energy must be conserved, the mechan- 
ical luminosity of the wind may not exceed the radiative 
luminosity of the star. The maximum mass loss rate that 
can be driven to infinity is therefore given by 
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where Mth- is the photon-tiring limit of the mass loss 
rate. If the mass loss rate at the stellar surface exceeds 



this limit, part of the mass will fall back toward the 
star. For a more thorough explanation of the analytical 
background of the porosity length formalism, as well as 
the continuum driving force, see OGS. 



Numerical approach 

We simulate the continuum driving by adding an ac- 
celeration term with the scaling of Eq.|2]to the momen- 
tum equation. The effect of photon tiring is included by 
reducing the local radiative luminosity in each gridcell 
by the amount of work required to drive the matter to 
that distance: 



(4) 



with L{r) the available local luminosity and 
M{r') — 4-nr'^p{r')v{r'). For a steady wind solu- 
tion, M(r') is a constant value, but if the star exceeds 
the photon-tiring limit this will no longer be the case. 
Note that Eq. |4] works two ways. If the local velocity 
is positive, the luminosity is reduced. If the velocity is 
negative, the integral term becomes negative as well and 
energy is added to the radiation field. In this equation 
we ignore for simplicity that the amount of work done 
by gas pressure is small compared to what is done by 
radiative driving. 



NUMERICAL RESULT 

For our sample simulations we use a 5QM,:;, SQRp) star 
with a surface temperature of 50 000 K. We ran several 
simulations for different luminosities to explore the be- 
haviour of continuum driven stellar winds. 



Below the photon tiring limit 

Figures [T] and |2] show the radial velocity and density 
of the wind for F — Lj^/LEdd = 3 and 6 respectively, both 
with and without the effect of photon-tiring included. As 
can be seen in these plots, the difference is comparatively 
small for the lower luminosity, reducing the wind veloc- 
ity by about 25%. However, for the high luminosity, the 
change is considerable. Here the terminal velocity of the 
wind is reduced by more than 50% and the wind velocity 
actually decreases somewhat in the outer layers as the lo- 
cal radiation field becomes too 'tired' to counterbalance 
the force of gravity. In both cases the terminal velocity 
of the wind is on the order of the escape velocity at the 
stellar surface. 

The mass loss rates for both winds is high. For F = 3 
the mass loss rate is about 0.04 Mq yr^K For F = 6 it 




FIGURE 1. Wind velocity (solid lines) for a T = 3, ce = 0.5 
simulation of a continuum driven wind from a SOMq, SORq 
star. The photon tiring effect only makes a small difference in 
the terminal velocity. 
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FIGURE 2. Similar to Fig. □ but with T = 6. For this lumi- 
nosity the photon tiring effect changes the wind velocity con- 
siderably. Note that the velocity actually decreases somewhat 
in the outer layers, as the local radiation force becomes less 
than gravity. 



is approximately 0.2 Mq yr^ . Such mass loss rates are 
several orders of magnitude more than can be achieved 
through line driving. Note that photon tiring does not 
change the mass loss rate. This depends exclusively on 
the porosity parameters and luminosity at the sonic point 
(OGS), which is defined as the point where grad = ^grav 
Only if this mass loss rate exceeds the photon-tiring 
limit does the mass loss rate driven to infinity depend 
on whether the photon tiring effect is included. 
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FIGURE 3. Similar stellar mass and radius as in previous 
figures, but luminosity at F = 10 in order to produce a surface 
mass loss rate of M = 2.34Mtii-. The line contours show the 
motion of individual mass elements, while the surface contour 
in the background shows the local radiative luminosity; both 
as functions of time and radius. The mass elements that leave 
the surface fall back down almost immediately due to a lack of 
radiative force. Only part of the mass progresses to larger radii, 
where a similar process repeats itself. Whenever an element 
moves outward, the radiative luminosity decreases causing ele- 
ments further out to slow down and/or fall back. If an element 
falls back, the radiative luminosity increases, which causes el- 
ements further out to accelerate. Note that the motion of mass 
elements reverses whenever radiative luminosity is particularly 
low or high. 



For these simulations the porosity length was set to the 
hydrodynamical scale-height at the stellar surface. 



Exceeding the photon tiring limit 

By increasing the luminosity to F = 10 and reducing 
the porosity length we increase the mass loss rate at the 
stellar surface to more than twice the photon tiring limit 
(Eq. 78 from OGS). The result can be seen in Fig. [3] 
which shows both the motion of individual mass ele- 
ments (lines) and radiative luminosity (surface contour) 
as a function of time. Mass is launched from the stellar 
surface, but most of it quickly falls back as the luminosity 
can not support it against gravity. Only a fraction of the 
mass that was initially launched manages to continue. At 
larger radii it encounters previously ejected mass, which 
is falling down again. Each time a mass element falls 
down, it re-energizes the radiation field (Eq. |4]l which 



allows the radiation to drive the mass element above fur- 
ther outward. If a mass element moves outward the radi- 
ation field loses energy, which causes the mass element 
above to decelerate or even fall back down. 

This process repeats itself several times, creating lay- 
ers of mass moving back and forth. Only part of the mass 
eventually reaches the point where its velocity is larger 
than the local escape velocity, allowing it to escape from 
the gravity of the star Clearly, the resulting wind has no 
stable solution. The mass loss rate of this wind is very 
high, reaching 80% of the photon tiring limit. However, 
the wind velocity is very low, since most of the energy is 
being used to lift this large amount of mass against the 
stellar gravity. 
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CONCLUSIONS AND PLANS FOR THE 
FUTURE 

We have shown that continuum driving can cause a star 
to lose large amounts of mass in a short period, provided 
that the star exceeds the effective Eddington limit. As 
long as the mass loss rate stays below the photon tiring 
limit, continuum driving can produce a steady wind with 
a high density and velocity on the order of the escape 
velocity. Once the mass loss rate exceeds the photon 
tiring limit the wind no longer has a steady solution. 
Instead the space close to the star becomes a tangle 
of mass elements moving up and down, depending on 
the amount of radiative energy that is available to drive 
them. The effective mass loss rate (that fully escapes the 
star) is still very large, but the wind velocity becomes 
quite low, since most of the available energy is used to 
lift the mass against gravity. Since continuum driving 
does not depend on the presence of metals it provides 
a very effective means to remove mass from zero- and 
low-metallicity stars. Such stars are also thought to be 
generally more massive than stars at solar metallicity 
(Schaerer [4], O'shea & Norman |2]), which makes them 
more likely to exceed the Eddington limit. 

In the future we intend to include 2D radiative hydro- 
dynamics in our simulations. We also want to compare 
our results with known super-Eddington objects in the 
nearby universe, such as LBVs and Novae. Hopefully 
this will allow us to create a framework for simulations 
of the mass loss rate of low metallicity stars. 
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